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Abstract 

The Pierre Auger Observatory has been designed to study the highest-energy 
cosmic rays in nature {E > 10^^'^ eV). The determination of their arrival direction, 
energy and composition is performed by the analysis of the atmospheric showers they 
produce. The Auger Surface Array will consist of 1600 water Cerenkov detectors 
placed in an equilateral triangular grid of 1.5 km spacing. The aim of this paper is 
to show that the addition of a "small" area of surface detectors at half or less the 
above mentioned spacing would allow a dramatic increase of the physical scope of 
this Observatory, reaching lower energies at which the transition from galactic to 
extragalactic sources is expected. 
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1 Introduction and Motivations 



Cosmic rays are observed in a wide range of energies spanning more than eleven decades, 
from energies below 1 GeV up to more than 10^° eV. At lower energies the spectrum follows a 
simple power law with an exponent equal to -2.7, compatible with galactic supernova remnant 
acceleration of charged nuclei. The spectrum slope becomes steeper at ~ 3 — 5 x 10^^ eV, a 
feature known as the knee, where the spectral index changes from -2.7 to -3.1. A spectrum 
compilation from 10^^ eV (just below the knee) to the highest detected energies is shown 
Fig. 1. This figure clearly depicts a large systematic uncertainty in energy calibration of 
~ 30% among the different experiments (and even larger for the Yakutsk data). Still, when 
the different experimental data are energy-renormalized there is fair agreement in shape up 
toE< 10^9-5 eV. 

The KASCADE experimental results [1,2.3] show evidence that the first knee is essentially 
caused by a decreasing flux of light primaries: as the atomic number increases, each element- 
associated knee moves to a higher energy. The KASCADE data in Fig. 1 show that the 
experiment is finely designed to inspect the knee region, i.e. in the ~ 10^^ — 10^^"^ eV energy 
range. Although it collects data above this energy and up to 10^^ eV, it does it with larger 
statistical uncertainties as it operates on the verge of its acceptance. 

In the region from ~ 10^^ ° to 10^^-^ eV the spectrum is reported to show two further traits: 
a break, called the second-knee, and a broad feature known as the ankle. The second-knee 
feature has been suggested to be a realization of the knee for the heaviest stable elements, i.e., 
Fe (see for example [4], [5]). It may represent the end of the efficiency of supernova remnant 
shock waves as accelerators or a change in the diffusion regime inside our Galaxy. Fig. 1 tells 
that AKENO [6], Yakutsk, with its latest analysis [7], and Fly's Eye stereo observe a clear 
second-knee-like feature. The figure also shows that KASCADE, Haverah Park, AGASA, 
HiRes and Auger have not enough acceptance at the second- knee region. HiRcs/MIA is at 
the limit of its acceptance, but there is still a hint of the second knee in its data. 

This acceptance gap will be filled by KASCADE-Grande, which will have an exposure area 
ten times larger than Kascade and will attain good statistics up to close to 10^^ eV by 2009, 
with composition studies up to (0.7 - 0.8) x 10^^ eV [8]. If so, it will bridge the undoubtedly 
very important energy range from the first to the second knee, and hopefully cast light on 
primary cosmic ray composition. 

As already mentioned, the second knee might be the end of the stable elements of the cos- 
mic ray Galactic component dominance [4,5,9]. Therefore, at higher energies an additional 
component would be required to account for the observed flux. In Refs. [10,11] it is conjec- 
tured that this additional component arises from extragalactic protons, together with some 
galactic heavy nuclei to conform a broad second knee feature. 

At energies beyond the second-knee feature, the spectrum clearly exhibits a broad depression, 
called the ankle. As shown in Fig. 1, the ankle has been observed around 3 x 10^® eV by 
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Fly's Eye [12,13,14], Haverah Park [15], HiRes [16], and it is not contradicted by Auger 
[17] at its preliminary stage. These results have been confirmed by Yakutsk [7] [18] and 
AGASA [19] but they locate it at a higher energy, around 10^^ eV. There are at least two 
physical interpretations of the ankle, intimately related to its nature: it may be the transition 
between the Galactic and cxtragalactic components [20,21,22] or the result of pair creation by 
extragalactic protons in the cosmic microwave background [10,23]. In the former model it is 
assumed that heavy nuclei are accelerated to the ankle energy within the Galaxy, whereas in 
the latter only up to the second knee and therefore the importance is not laid on the ankle but 
on the broad trait which is interpreted as a pair-production dip. It has also been suggested 
that the ankle is the result of diffusive propagation of extragalactic nuclei through cosmic 
magnetic fields [24,25]. To pinpoint the correct model both a reliable, high quality spectrum 
with a well calibrated absolute energy and a detailed composition study are required, as 
the pair-production dip model relies on cosmic rays being essentially protons beyond the 
second-knee [11,22]. 

In any case, the energy region spanning from log{E/eV) ~ 17.5 to ~ 19.0 very likely com- 
prises the transition from Galactic to extragalactic cosmic rays. 

1.1 The Pierre Auger Observatory 

As of today the forefront experiment in the ultra-high energy cosmic ray arena is the Pierre 
Auger Observatory [26], that aims at building two observatories, one in each hemisphere. 
The construction of the Southern Observatory started in 2000. It is already taking data 
and approaching completion. It is located in the region called "Pampa Amarilla", close to 
Malargiie, at the south of the Province of Mendoza in Argentina (35.0° — 35.3°5', 69.0° — 
69.3°VF), therefore having a full view of the Galactic center and its surroundings. Auger's 
two distinctive features are its exceptional size and its hybrid nature. It spans over an area 
of 3000 km^ and is constituted by a surface array of 1600 water Cerenkov detectors placed 
on a 1.5 km triangular grid plus 24 fluorescence detector telescopes placed in 4 buildings on 
the Surface Detector (SD) array periphery and overlooking it. Consequently, it will provide 
a large number of events with low systematic uncertainties. 

Pertaining to this work is the Auger surface detector array. It has a versatile trigger system 
[27] designed to operate in a wide range of primary energies and arrival directions. Apart of 
the customary Signal-over-Threshold (ST) trigger. Auger bases its SD local trigger system 
on a Time-over-Threshold (ToT) trigger, which requires the signal to be above the ToT 
threshold (which is much lower than the ST threshold) during at least 325 ns in a 3 jis 
interval. The array trigger requirement is satisfied when at least 3 surface detectors in a 
compact configuration detect a local trigger in coincidence. With this set of conditions the 
Southern Observatory is fully efficient above 3 x 10^^ eV [28]. For energies below this value, 
the detection efficiency decreases rapidly and is composition dependent. Also, the shower's 
parameters (i.e. energy, zenithal angle, and core position) reconstruction is degraded since an 
average of three tanks does not suffice to convincingly sample the shower lateral distribution 
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function (LDF). Detection efficiency also becomes fluctuating with atmospheric conditions 
(temperature and pressure), introducing unknown uncertainties in sky coverage. 

1.2 The Graded Infill 

In this paper it is proposed to enhance the Auger acceptance down to 10^^ eV by means of 
a graded infill of surface detectors deployed at smaller distances, over an area much smaller 
than Auger because of the much higher fiux at lower energies. This infill array would use 
the same technology of Auger surface detectors and benefit from the existing knowledge 
and infrastructure at the site (detector design, data acquisition, analysis tools, etc.). Such 
an extension would allow: (a) to make a reliable measurement with low reconstruction un- 
certainty of the energy spectrum from the second knee to the ankle, (b) to make more 
accurate anisotropy measurements at lower energies avoiding the many spurious effects like 
the temperature dependence of the event rate and the composition dependent efficiency, and 
(c) to bridge with a sizeable overlap the data between KASCADE-Grande and Auger, thus 
increasing the reliability of the results obtained. 

An infill would also give an experimental handle on fluctuations at the verge of Auger 
acceptance, i.e. ~ 3 x 10^* eV, as it will allow event reconstructions with two or more 
non-overlapping subsets of detectors at 1500 m distance. In a similar fashion it would permit 
a very precise experimental determination of the acceptance of the full array, which will also 
serve to validate detector simulation packages. Additionally, the inflll could contribute to a 
better primary discrimination through a comparison with the dependence of the acceptance 
on composition. 

Three different infill configurations are considered, which can be obtained by adding surface 
detectors to the Auger grid, which has a spacing of 1500 m. These infills give new triangular 
arrays with spacings between neighbouring detectors of, respectively, 866, 750, and 433 m. 
In order to assess the impact of such an enhancement, we study by means of simulations the 
dependence of the resolution for different shower parameters (arrival direction, core position, 
lateral distribution of shower particles, energy, etc.) on detector spacing, for different primary 
energies, composition, and zenith angles. We will conclude that in order to reach the above 
mentioned energy range with unitary efficiency, a graded infill is required, i.e, an array with 
a 750 m-grid and, in an even smaller area, an infill of this infill, with a spacing of 433 m 
between SDs. 

This paper is organized as follows: section 2 describes the simulations of showers and detector 
response, as well as the analysis methods used. Section 3 contains our acceptance calcula- 
tion for different infill configurations. In sections 4 to 6 the improvement on reconstruction 
of shower parameters is analyzed. Some considerations about the possible contribution to 
composition studies are given in section 7; conclusions are presented in section 8. 
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2 Simulations 



For this work we generated a library of extensive air showers using the Monte Carlo simulation 
code Aires 2.6.0 [29] (see also [30]). Two types of primaries (proton and iron) were considered, 
arriving with three characteristic zenith angles: 0°, 30° and 45°, with five different energies: 
10^^-5 eV, 10^^-^^ eV, 10^^ eV, lO^^-^^ eV and 10^^-^ eV. For each energy, zenith angle and 
primary composition, 100 showers were simulated with a uniform azimuthal distribution. 
With each of these showers, 5 events were generated by triggering the surface array at 
random impact points. 

Using the code SDSim (v3r0) [31], the response of an Auger-like surface array was simulated, 
consisting of 37 detectors covering a hexagonal area of 52 km^ (see Fig. 2, top left). Three 
different infill configurations were added to this 1.5 km triangular grid of detectors: 

• a) One detector at the center of each triangle of 1500 m, resulting in a triangular grid with 
a spacing of 866 m (Fig. 2, top right) 

• b) Detectors at half the distance between Auger detectors, resulting in a triangular grid 
with a spacing of 750 m (Fig. 2, bottom left) 

• c) Detectors at the center of each triangle of 750 m from grid b), resulting in a triangular 
grid with a spacing of 433 m (Fig. 2, bottom right) 

All Auger trigger levels [27] with their relevant signal thresholds have been used in these sim- 
ulations. After all trigger requirements are met, a full reconstruction of the arrival direction 
and the lateral distribution function is performed using the event reconstruction package 
Er(v3r4) provided by the Central Data Acquisition Group of the Auger Observatory [32]. 
Reconstructions invoke a Lateral Distribution Function (LDF) parameterization which de- 
scribes the expected integrated signal at a distance r from the shower core. We chose the 
functional form: 



where Tq =1000 m and Siqqo, (3 and 7 are parameters to be adjusted. 

At the lowest energy considered (3 x 10^^ eV) the number of showers that are able to trigger 
the Auger 1500 m-array is marginal and, correspondingly, the determination of the resolution 
of shower parameters is physically meaningless. 

The array acceptance and the arrival direction and LDF reconstruction accuracy depend 
strongly on the number of stations with signal above threshold. Fig. 3 shows the number 
of triggered stations versus primary energy for different infill spacings. It is seen that, as 
the shower area on the ground increases with increasing energy, the number of triggered 
detectors gets correspondingly larger. For a given detector spacing, this number is expected 
to be proportional to E"^/^, where (3 is the slope of the LDF. The values and fits shown in 
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Fig. 3 correspond to iron primaries. 



3 Acceptance Determination 

The shower detection efficiency of a surface array depends crucially on the spacing between 
detectors and on the trigger configuration, i.e., the number of detectors that are required 
to have signal above a certain threshold for the event to be recorded. The surface detector 
array of the Pierre Auger Observatory, with a spacing of 1500 meters between neighbouring 
detectors, is nearly fully efficient for showers of more than 3 EeV, if a 3-detector trigger 
configuration is required [27]. This value has been obtained by Monte Carlo simulations of 
the detector response. 

The instantaneous acceptance of the surface array is given by: 

Acc{E)n+ = J PN+{f, E, 6) ■ cos{e) ■dr^-27r- sin(^) • dO (2) 

where Pjv+ is the probability of having a positive trigger with at least N detectors for a 
shower of energy E hitting the ground in the position f, with zenithal arrival direction 9. In 
the full acceptance regime, i.e. when P/v+ = 1, the acceptance of an array of area A, and for 
a maximum zenith angle of 9max, is Acci = At: sm^{9max)- For what follows, 9max is set to 
60°. 

To calculate the probability Pn+, we need to consider the local trigger probability for a single 
detector (say, the i-th detector of the array). To a good approximation, this probability is 
determined solely by the total expected integrated signal of the station. Thus, if the i-th 
local station is located at fj and has an expected signal Si, its trigger probability is given by 

P{Si)^P{LDF{\fi-r\,E,9)). (3) 

Here we are making the general assumption that the expected signal at a certain distance 
from the shower core r can be modelled by a Lateral Distribution Function (LDF) (see 
section 1). 

To obtain P{S), a set of showers in the 10^^ - 10^^ eV energy range with zenith angles up to 
60° were simulated, together with the corresponding detector response to obtain the signals 
in the detectors. These simulated signals, which include the fluctuations both in the shower 

development as well as in the detector response, are used to reconstruct the shower core 
position, energy and incidence angle by adjusting an LDF to the signals in the triggered 
detectors. 

From a large number of simulated showers with different core positions the P{S) can be 
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inferred as 



where Nt+{S) is the number of triggered stations with an expected signal S, and Non{S) is 
the total number of active stations in that signal bin. Given this P{S), the integral (2) can 
be calculated numerically. 

Fig. 4 shows the resulting relative acceptance Acc/Acci as a function of energy, for pro- 
tons {open symbols) and iron {dosed symbols) and for the different detector spacings under 
consideration. It can be observed that a separation of 750 m between detectors guarantees 
a detection efficiency of 95% for proton (iron) showers with an energy of 3.6 x 10^^ eV 
(1.8 X 10^'^ eV). With a detector spacing of 433 m, a similar acceptance can be obtained for 
showers of as low as 9 x 10^^ eV and 4 x 10^^ eV for proton and iron primaries, respectively. 
According to this calculation, at 3 EeV the Pierre Auger array would attain full acceptance 
for a 3-detector trigger, irrespective of composition, in agreement with ref. [28]. 

Considering a differential flux of cosmic rays ^ following a power law with spectral index 
—2.84 as quoted by Auger in [17]: 

dd) f E \ -^-^ ,911 



the number of expected events in one year, for events with energy larger than Eq with a 
zenith angle below 9 max — 60°, for a detector of area A is given by: 



= 795 ■ 



■ A - 
.20km?. 



t 

yr 



.lEeV. 



-1.84 



(6) 



Consequently, an infill with an effective area of ~20 km"^ (for well-contained events), which 
would comprise a total of 64 surface detectors separated by 750 m, could accumulate a 
statistically significant number of events with a relatively small effort: only 44 detectors 
beyond the regular Auger array would be required, which amounts to less than 3% of the 
full Auger Observatory. These detectors would operate in an energy region in which Auger 
ceases to be fully efficient (which renders the reconstruction not reliable), down to energies 
in which the regular Auger array does not trigger any more. 

As already mentioned, to detect showers with full efficiency at energies as low as 10^^ eV, a 
denser infill, with a detector spacing of 433 m, is required. A statistically significant number 
of events can be recorded with the addition of only 10 detectors with such a spacing, covering 
an effective area of 1.6 km^. 
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4 Angular Resolution 



The accuracy in the determination of the arrival directions of cosmic rays is fundamental in 
the search for their origin and in the study of anisotropics, although it is noted that only 
neutral particles are not deflected during their propagation in magnetic fields. We define the 
angular reconstruction uncertainty as the space angle © subtended by the real {Ureal) and 
reconstructed (Rrec) directions, being cos{Q) — Rreai • Rrec- 

In Fig. 5 we present the 68% confidence level for the arrival direction reconstruction uncer- 
tainty, AG, as a function of detector spacing and shower energy. Results are presented for 
iron {left) and proton {right) at different injected energies. It can be seen that the angular 
resolution for events detected by a 750 m infill array will be 2 times better than for a 1500 
m-array, giving values close to 1° for proton and iron at 10^^ eV. 

It is worth analyzing the dependence of A0 on the number of triggered detectors N^et, as 
shown in Fig. 6 for proton and iron primaries and for all considered values of primary energy 
and detector spacing. As expected, as the average number of triggered stations participating 
in the reconstruction increases, the resolution of the arrival direction improves. This behavior 
is seen for both primaries considered. The full lines represent a fit to the data for proton 
{blue), iron {red) and for both {green). The best fit for the complete set of data is given by 
Ae = 5.8° X + 0.6°. 

As the trigger conditions used in our simulations are the same as for the Auger surface 
detector [27], the dependence of the angular resolution on the number of triggered stations 
for the 1500 m-array may be compared to that obtained by the Auger SD [35]: the values 
of A© corresponding to 3, 4 and 5 stations (2°, 1.6° and 1.2°, respectively) are in good 
agreement with the results presented in [35], where it is found that for the 3-fold events. 
Auger has an angular resolution of about 2°, for 4- fold events of 1.7° and for 5-fold events 
of 0.9°. 

It is worth comparing the obtained angular resolution with that quoted by other experiments 
designed to operate in a similar energy range. The 750 m-infiU is expected to give an arrival 
direction resolution of 1° at 10^^ eV while AKENO achieves a pointing accuracy of 3° [36]. 
At 10^'' eV, Kascade-Grande [4] is expected to have an angular resolution of 0.3° and IceTop 
[37] claims an accuracy of 0.6°. These values should be compared to the proposed 433 m-infiU 
which will operate in this energy region and will have an angular resolution smaller than 
0.7°. 
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5 Core Position Resolution 



A very important parameter for the reconstruction of shower geometry and energy is the core 
position, i.e. the intersection point of the shower axis with the ground. A first estimation is 
given by the barycenter of the 3 highest-signal detectors. This value is used as input for an 
iterative fitting process in which the lateral distribution of the detector's signal is adjusted. 

To estimate the expected uncertainty on core reconstruction for different configurations, we 
evaluate the difference between the real and reconstructed core positions: 



where Xrec, Urec are the coordinates of the shower core (on the array plane) obtained from 
the full reconstruction process, and Xreah Ureal are the input values for each of the simulated 
showers. We define the core position resolution (for each energy, angle and primary type 
considered) as the value at which the integral of the AC ore distribution is 68% of the total 
integral. 

The main results of this section are shown in Fig. 7. Since the lateral distribution function 
decreases rapidly with distance, the presence of stations near the core plays a fundamental 
role in its position determination. This situation is favored with smaller spacings between 
detectors: as the distance between detectors decreases, the core position is better determined. 
Note that for the Auger spacing of 1500 m the empirical RMS value found in [33] is ~102 
m, which is in good agreement with our results. Shower development for proton primaries 
present more fiuctuations than iron and this is refiected producing bigger errors. Almost the 
same behavior as depicted in Fig. 7 for showers with 9 — 30° was observed for vertical and 
more inclined (45°) showers. 

At 433 m spacing, the number of triggered detectors is large enough (~ 20, see Fig. 3) to 
have a core resolution of ~ 15 m irrespective of the considered energy and primary type. 
Regarding the 750 m infill, it is seen that the core position resolution improves by a factor 
of ~ 4 with respect to Auger. This will have an important impact on the reconstruction, 
giving a better angular and energy resolution, as the LDF fitting is also more accurate. 

Fig. 8 shows the core position resolution as function of the average number of triggered 
detectors, for both primary types, and all considered energies and spacings, at ^ = 30°. 
The lines represent the best fit to the values obtained for the different primaries {blue for 
proton and red for iron) and for both (green). The core position dependence with the average 
number of triggered stations can be adjusted by ACore = 820 m xN^Jj.'*'^ + 13 m. 




(7) 
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6 S(600) Determination 



For the events registered by a ground array of surface detectors, the determination of the 
energy consists of two steps. First, an energy dependent parameter called S'(ro) is assigned to 
a given shower. This parameter is the time-integrated signal expected in a surface detector 
placed at a distance Tq from the shower axis. Its value is obtained by fitting an empirical 
LDF, with a predefined functional form, to the observed lateral signal distribution and 
interpolating its value at the distance Tq from the core. The second step involves an energy 
calibration of 5'(ro). In a hybrid detector such as Auger this can be performed either by 
resorting to hybrid events, for which both the fluorescence and surface data can be well 
reconstructed independently [32] , or by obtaining the conversion function from Monte Carlo 
simulations. 

The choice of the parameter Tq for the energy conversion is directly related to detector 
spacing and primary energy considered. It has been proven [33] that Tq = 1000 m is the 
optimal value for a grid with 1.5 km spacing and for energies larger than ~ 5 x 10^^ eV. 
However, since our goal is to reach lower energies with a smaller detector spacing, the chosen 
parameter is instead 5'(600), i.e., the time- integrated signal expected in a detector at 600 
meters from the shower core. 

The accuracy in the determination of the parameter S{ro) depends on the detector resolution 
and on sampling fluctuations. Additionally, this accuracy has also a strong dependence on 
shower to shower fluctuations which are intrinsic to the shower development and can not be 
eliminated, as well as on shower direction and relative core position [33]. Moreover, the non- 
linear interpolation used in the determination of 5'(ro) increases even further its uncertainty. 
It was found empirically [34] that for the Auger grid the statistical uncertainty on S'(IOOO) 
reaches 10% (corresponding to an error of 50 m for the core location), the same order of 
magnitude as the sampling contribution. These values were obtained for showers with E 
> 5 X 10^^ eV. 

In Fig. 9 we plot the ratio a5'(6oo)/>S'(600), where (Ts(mo) is the dispersion of the reconstructed 
S(600) distribution, for both primaries: iron (triangles) and proton (circles), for different 
injected primary energies. When the distance between detectors is reduced by a factor of 
two, the accuracy of the new energy parameter S(600) for E — 10^^ eV is ~ 10% for iron 
primaries whereas for protons this value is about 13 %. 



7 Primciry Composition 

In the same way as the cosmic ray's energy and arrival direction are better determined us- 
ing the information provided by an inflU array, it is expected that other variables carrying 
information about the identity of the primary particles would behave similarly. For example. 
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for surface detectors, it is known that the number of muons at a certain distance of the 
shower core, the slope of the LDF, the arrival time profile of shower particles, and the radius 
of curvature of the shower front can be employed as indicators of the primary composition 
[38] . The infill allows not only a more precise shower reconstruction as shown in the previous 
sections, but also helps towards composition analyses down by another decade in energy, 
completely covering the second knee and the ankle. As an example, we studied the distri- 
bution of the LDF slope and the shower front curvature for different characteristic energies 
and both primary types, with shower and detector simulations for 9 — 30°. 

An estimator of the power to discriminate different primaries is given by the parameter 77 
defined as: 

^ _ \Med{p) - Med{Fe)\ 



where Med{p) and Med{Fe) are the medians of the distributions for proton and iron respec- 
tively, apfFe^^^ the dispersions of the distributions [38]. 

Fig. 10 (top) shows the distribution of the slope /? of the LDF function (see Eq. (1)), for 
both primaries and for two different spacings: 1500 m {left) and 750 m {right), at a primary 
energy of 3.0 EeV, i.e. at the lower limit of Auger full trigger efficiency. An inspection of 
this top panel shows that /3 is quite inappropriate for primary discrimination for the current 
Auger spacing but it seems useful for the 750 m infill array. 

Fig. 10 (bottom) shows the P distributions for 10^^ eV showers on the Auger surface array 
{left) and 10^^ eV showers on the 750 m infill array {right). This bottom panel shows that 
the 750 m infill has a similar discrimination power as the full Auger array, but at an energy 
one order of magnitude lower. 

Performing the same kind of analysis with the radius of curvature of the shower front, R, 
yields similar results. In Fig. 11 (top) the distributions for 3 x 10^^ eV are presented. Once 
again, at this energy the infill array helps to statistically discriminate between proton and 
iron, whereas this is seriously compromised for the Auger spacing. The bottom panel shows 
that the power of discrimination at 10^^ eV using the infill information is slightly better than 
that of the present Auger spacing for 10^^ eV. 

The improvement in these parameters gives a clear hint that composition studies would 
greatly benefit from the information provided by an infill array. Moreover, a multiparametric 
analysis as suggested in [38], which weights in a large number of variables, would allow a 
much more precise composition resolution. 
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8 Conclusions 



In this work we showed that enhancing the Pierre Auger Observatory by increasing the 
density of surface detectors in a small area (20 km^ of the total 3000 km^) would bring 
considerable advantages to the study of high energy cosmic ray physics, at a relatively 
low cost (3% of the complete Auger SD). From the operational point of view, such an 
enhancement could be smoothly integrated to the existing surface array. 

With an infill array much could be gained in shower reconstruction at the lower limit of 
Auger energies: a better LDF fitting due to a larger number of triggered detectors, a 2x and 
4x improvement in arrival direction and core position resolution respectively and a relative 
uncertainty of less than 15% in 5(600) determination. 

Also, an infilled array with 750 m spacing is equivalent to 3 superimposed Auger-like arrays of 
1500 m spacing. Single events can then be reconstructed with 3 different subsets of detectors 
at Auger spacings, allowing a direct study of fluctuations and uncertainties. The data coming 
from an infill array would thus serve to check the behavior of the original array and to validate 
further the end-to-end simulation and reconstruction processes. 

The extension of the detector's full efficiency interval down to ~ 10^^ eV to englobe com- 
pletely the second knee and the ankle regions would be difficult to overestimate, allowing 
to test in an unprecedented way competing models of Galactic and extragalactic cosmic 
ray production and propagation. To this end, better composition analyses could be per- 
formed since distinctions among theoretical models largely depend on the primary type. 
Furthermore, an additional independent cross-check could be achieved by overlapping with 
KASCADE-Grande in the lower energy region. 
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Fig. 1. Differential cosmic ray spectra multiplied by E^, showing the knee, the suggested second 
knee and the ankle features. Compilation from KASCADE [3], AKENO [6] Yakutsk [7], Haverah 
Park [15], Fly's Eye Stereo [39], HiRes-MIA [40], HiRes I and II [40], AGASA [19], and AUGER 
[17]. The KASCADE data was obtained using the QGSJET 01 interaction [41] and AKENO data 
is from the 1 kra^ array. 
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Fig. 2. Different possible infill configurations generated inside the original array with 1500 m spacing 
{top left): a detector placed at the center of each triangle {top right) makes a grid of 866 m. A 
detector added at half distance between the original detectors (bottom left) results in a grid of 750 
m. A detector placed at the center of each triangle of 750 m (bottom right) produces a grid of 433 
m. 
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Fig. 4. Relative acceptance for different infill configurations for proton (open symbols) and iron 
(closed symbols). 
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Fig. 5. Space angle uncertainty as a function of detector spacing for iron (left) and proton (right) 
for 6 = 30°. Lines are drawn only to guide the eye. 
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Fig. 6. Space angle uncertainty as function of the mean number of triggered stations for both 
primaries: proton (circles) and iron (triangles) for all considered energies and detector spacings. 
These results correspond to ^ = 30°. 
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Fig. 7. Core position resolution as a function of detector spacing for both primaries: iron (left) and 
proton (right) for 6 = 30°. Lines are drawn only to guide the eye. 
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Fig. 8. Core position resolution as a function of the mean number of triggered stations for both 
primaries: proton( circ/es) and iron {triangles), for aU considered energies and detector spacings and 
9 = 30°. 
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Fig. 9. Relative dispersion on S(600) defined as 0-5(600) /5'(600) for iron (triangles) and proton 
(circles) as function of input energy for 750 m spacing and for 6 = 30° 



19 



20 - 3EeV- 1.5 km 
11=0.1 



i£L 



-5 -4.5 -4 -3.5 -3 -2.5 -2 
P 




3 EeV - 0.75 km 



ri= 0.9 




-5 -4.5 -4 -3.5 -3 -2.5 -2 
P 



1 EeV - 0.75 km 



11=0.5 



Proton ■ 
Iron - 



-5 -4.5 -4 -3.5 -3 -2.5 -2 



Fig. 10. Top: Distributions of /3 (slope of the LDF) for iron and proton at the hmit of Auger 
acceptance, i.e. 3 x 10^*^ eV, for a spacing of 1500 m (left) and an infill of 750 m (right). Bottom: 
Distributions of (3 for both primaries at 10^^ eV for a Auger 1500 m-array {left) and at 10^^ eV for 
a detector spacing of 750 m (right). 
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Fig. 11. Top: Distributions of the curvature radius of the shower front for iron and proton at the 
hmit of Auger acceptance, i.e. 3 x 10^^ eV for Auger spacing [left] and an infih of 750 m (right). 
Bottom: Distributions of the curvature radius of the shower front for both primaries at 10^^ eV 
for a Auger 1500 m-array (left) and at 10^*^ eV for a detector spacing of 750 m (right). 
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